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We develop a relativistic mean field (RMF) description of deformed nuclei with pairing 
correlations in the BCS approximation. The treatment of the pairing correlations for nuclei 
whose Fermi surfaces are close to the threshold of unbound states needs special attention. 
With this in mind, we use a delta function interaction for the pairing interaction to pick up 
those states whose wave functions are concentrated in the nuclear region and employ the 
standard BCS approximation for the single-particle states obtained from the RMF theory 
with deformation. We apply the RMF + BCS method to the Zr isotopes and obtain a good 
description of the binding energies and the nuclear radii of nuclei from the proton drip line 
' to the neutron drip line. 

O 

■ 81. Introduction 

It is our strong desire to obtain a model valid for all nuclei, including unstable 
ones, from the proton drip line to the neutron drip line. The relativistic mean field 
(RMF) theory has been used to describe such nuclei with one parameter set in all 
mass regions yy e h ave t Q include deformation and pairing correlations into the 
RMF model for a suitable description of finite nuclei. There is an extended study 
of all the even-even nuclei over the entire mass region by Hirata et al. including 
only deformation.® This calculation provides a good account of all the nuclei and 
indicates that almost all nuclei, except for those with most of the magic numbers, are 
deformed. Their calculation, however, does not include pairing correlations, because 
the conventional BCS treatment with the constant pairing interaction is not able to 
treat the case in which the Fermi surface is close to the unbound threshold.® 

The pairing correlations, on the other hand, have been treated nicely in the 
framework of the relativistic Hartree-Bogoliubov (RHB) method by Meng et al.® 1 ® 
They were able to treat nuclei near the unbound threshold in the RHB framework. 
They solved the RHB equation in the coordinate space with box boundary condi- 
tions. This method has been applied to many proton magic nuclei and is capable 
of describing many interesting features, such as the giant halo where the neutron 
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density distribution extends far from the nuclear region. This method is, however, 
limited to spherical nuclei. When extended to deformed system, this method requires 
such great compuatation time that we are not able to make calculations for all the 
nuclei in the periodic table in a systematic manner. 

Recently, there was an interesting suggestion made by Yadav et al. that us- 
ing the delta function interaction in the BCS formalism with proper box boundary 
conditions can provide a good description of the proton magic nuclei.^ Indeed, the 
calculated results obtained from the RMF+BCS method and those obtained in the 
RHB framework are nearly equal for those proton magic nuclei with spherical shapes. 
The justification of this method was provided by Sandulescu et al. by solving the res- 
onance states and taking into account the width effect exactly in the non-relativistic 
Skyrme Hartree-Fock framework.^ The applicability and justification of such a delta 
function interaction is also discussed extensively in a work by Dobaczewski et alP 
and references therein, where it is shown in the context of HFB calculations that the 
use of a delta force in finite space simulates the effect of finite range interaction in 
a phenomenological manner and can take into account the effects of unbound states 
properly. It is very interesting, therefore, to apply this prescription to deformed 
nuclei. 

The self-consistent RMF theory was first extended to treat deformed nuclei by 
Price et alff^ and Gambhir et alP^ In Gambhir's work, the nucleon wave functions 
and the meson fields are expanded in terms of the harmonic oscillator wave functions. 
We employ this method for the mean field part and replace the BCS part with a 
constant pairing interaction by one with a delta function interaction. The fact that 
we use the nuclear wave functions (in particular, we compute the overlap between 
occupied and unoccupied states), allows us to pick up states in the continuum whose 
wave functions are concentrated in the nuclear region. This method, at least, is 
effective for spherical nuclei and may also be effective for the deformed case. We 
mention that the relativistic Hartree-Bogoliubov theory was worked out using the 
expansion method for deformed nuclei.^ 

In this paper, we formulate the RMF theory with the BCS method for pairing 
correlations. In §2, we present the RMF formalism with deformation and pairing. 
In §3, we apply the method to Zr isotopes from the proton drip line to the neutron 
drip line. We compare the calculated results with those obtained with the assump- 
tion of spherical shapes in RCHB and those obtained without including the pairing 
correlations. In §4, we provide the summary of the present work. 

§2. RMF with deformation and pairing 

We present here the formulation of the RMF theory with deformation and pairing 
correlations. We employ the model Lagrangian density with nonlinear terms for both 
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a and oj mesons, as described in detail in Ref., which is given by 



C = ipii^dn - M)ip + ^ad^a - ^m 2 a a 2 - ^g 2 (J A - \gz<J A - g^ipaip 
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(2-1) 

where the field tensors of the vector mesons and of the electromagnetic field take the 
following forms: 




d li f^-d u p^-2g p e abc P y v 



(2-2) 



8fj,A u 8 u A p , 



and other symbols have their usual meanings. 

The classical variational principle leads to the Dirac equation, 

[-iaV + V(r) + (3 (M + S(r))\ ^ = 

for the nucleon spinors and the Klein- Gordon equation, 



(2-3) 



{-A + m 2 a }a(r) 
{-A + ml}ojV( r ) 
{-A + m 2 p }p a ^(Y) 



-g a p s (r) - g 2 (J 2 (r) - # 3 o" 3 (r) 
^i M (r) +54^(r)a;^(r), 

ejp(r), 



(2-4) 



-AA»(r) 

for the mesons. Here, V(r) represents the vector potential 

V(r) = ^7 M ^(r) + ff^W(r) + e^V^W, 



and 5(r) is the scalar potential 



S(r) 



(2-5) 



(2-6) 



For the mean field, the nucleon spinors provide the corresponding source terms: 



Ps{r) 
J"(r) 
J aM (r) 

#(r) 
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E 
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(2-7) 



Here, the summations are taken over the valence nucleons only. It should be noted 
that as usual, the present approach ignores the contribution of negative energy states 
(i.e. no-sea approximation), which implies that the vacuum is not polarized. The 
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coupled equations (|2-3|) and (|2-4|) are non-linear quantum field equations, and their 
exact solutions are very complicated. For this reason, the mean field approximation 
is generally used; i.e., the meson field operators in Eq. (|2-3|) are replaced by their 
expectation values. In this treatment, the nucleons are considered to move indepen- 
dently in the classical meson fields. The coupled equations are solved self-consistently 
by iteration. 

The symmetries of the system simplify the calculations considerably. In all the 
systems considered in this work, there exists time reversal symmetry, so there are 
no currents in the nucleus and therefore the spatial vector components of uj^, p a ^ 
and A^ 1 vanish. This leaves only the time-like components, u>°, p a0 and ^4°. Charge 
conservation guarantees that only the 3-component of the isovector p 00 survives. 

2.1. Axially symmetric case 

The RMF theory was extended to treat deformed nuclei with axially symmetric 
shapes by Gambhir et alP^ To make clear the notation used, here a brief review of 
the RMF method for axially deformed nuclei is given. 

Many deformed nuclei can be described with axially symmetric shapes. In this 
case, rotational symmetry is lost, and therefore the total angular momentum, j, is 
no longer a good quantum number. However, the densities are still invariant with 
respect to rotation about the symmetry axis, which is assumed to be the z-axis in the 
following. It is then useful to work with cylindrical coordinates: x = r±_ cos <p,y = 
r±simp and z. For such nuclei, the Dirac equation can be reduced to a coupled 
set of partial differential equations in the two variables z and r±. In particular, the 
spinor tpi with the index i is now characterized by the quantum numbers l7j,7Tj and 
ti, where J?j = + m Si is the eigenvalue of the symmetry operator J z , TTi is the 
parity and t, is the isospin. The spinor can be written in the form 



fi(r) _ 
»fli( r ) J V2tt 



( /+(z,r ± )e i (^" 1 / 2 )^ \ 
fr(z,r ± )e i ^+ 1 / 2 ^ 
ig+iz^^e^- 1 / 2 ^ 
\ igi(z,r 1 _)e i ^+ 1 / 2 ^ ) 



XtM- (2-8) 



The four components f i (z,r±) and g i (z, r±) obey the coupled Dirac equations. 
For each solution with positive f2i, ipi, we have the time-reversed solution with the 
same energy, ife = Tipi, with the time reversal operator T = —ia y K (K being the 
complex conjugation). For nuclei with time reversal symmetry, the contributions to 
the densities of the two time reversed states, i and i, are identical. Therefore, we 
find the densities 

PsiV = 2Y J ((\ft\ 2 + i/n 2 ) t (is+i 2 + bn 2 )) (2-9) 

and, in a similar way, p^ and p c . The sum here runs over only states with positive 
These densities serve as sources for the fields (j) = cr, w°,p 00 and ^4°, which are 

determined by the Klein-Gordon equation in cylindrical coordinates. 

To solve the RMF equations, the basis expansion method is used. We closely 

follow the details, presentation and notation of RefP^ For the axially symmetric 
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case, the spinors and gf in Eq. (|2-8[) are expanded in terms of the eigenfunctions 
of a deformed axially symmetric oscillator potential, 



V osc (z,r ± ) = ^Muj 2 z z 2 + ^Mojlr\. 



(2-10) 



Then, imposing volume conservation, the two oscillator frequencies lo± and lo z can 
be expressed in terms of a deformation parameter, u> z = uiq exp ( — J j-fio ) and 



u± = u; exp g|/3o^ 

The basis is now determined by the two constants loq and (3$. The eigenfunctions 
of the deformed harmonic oscillator potential are characterized by the quantum 
numbers, \a) = \n z ,n r ,mi,m s ) , where mi and m s are the components of the orbital 
angular momentum and of the spin along the symmetry axis. The eigenvalue of J z , 
which is a conserved quantity in these calculations, is fi = mi + m s . The parity is 
given by ir = (_) n ^+ m i. 

The eigenfunctions of the deformed harmonic oscillator can be written explicitly 

as 

(2-11) 



with 



$ a (z, ri, <p, s, t) = <j> Uz (z)<f>% ( r± )^=e im ^xm s (s) X t a (t), 

V 2ir 



N, 



Arm 

n Lr? m i /2 jL rn i(7?)e -r,/2 ) 



(2-12) 



where £ = z /b z and r] = r\jb\. The polynomials H n {C) and L™(j]) are the Hermite 
polynomials and the associated Laguerre polynomials, as defined in Refp3 ;l The 
quantities N Uz and N™ 1 are normalization constants. 

The spinors ff 1 and g- in Eq. (|2-8jl are explicitly given by the following relations: 



/+(^i) 

/r( z > r -L) 



£ fa <PnA z )4>n r 
a 

a 

£ % Vn 

£ 

/3 



(ri.), 



(*)*£- 1/2) (n_), 



(2-13) 



The quantum numbers a max and /3 max are chosen in such a way that the correspond- 
ing major quantum numbers N = n z + 2n p + m; are not larger than iVir + 1 for the 
expansion of the small components and not larger than Np for the expansion of the 
large components. 

2.2. Pairing with delta function interaction 

Based on the single-particle spectrum calculated with the RMF method de- 
scribed above, we carry out a state-dependent BCS calculation.^'^ The gap equa- 
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tion has a standard form for all the single particle states, 

A k = -- V , %k ' Ak ' (2-14) 

where ey is the single-particle energy and A is the Fermi energy. The particle number 

condition is given by 2 ^ v\ = N. In the present work, we use a delta force for the 

fc>o 

pairing interaction, 

V = -Vo8(n - r 2 ), (2-15) 

with the same strength Vq for both protons and neutrons. The pairing matrix element 
for the 5-function force is given by 

% = (ii\V\jj) - (ii\V\jj) = -V J d 3 r [44^3 ~ > ( 246 ) 

with the pairing energy defined by impair = — Yl ^kUkVk- Equations l|2-3j) and (|2-4|) . 

the gap equations (|2-14jl . and the total particle number condition N for a given 
nucleus are solved self-consistently by iteration. 



§3. Numerical calculation for Zr isotopes 



We apply the formalism to the Zr isotopes from the proton drip line to the 
neutron drip line. For the RMF Lagrangian, we use the ^4-dependent parameter 
set TMAP^'C^'^ The parameter values are as follows. The masses of nucleon, 
<r, to and p mesons are, respectively, M = 938.900 MeV, m a = 519.151 MeV, 

= 781.950 MeV, m p = 768.100 MeV. The effective strengths of the couplings 
between various mesons and nucleons have the values g a = 10.055 + 3.050/^4 , 
g w = 12.842 + 3.191/A 04 and g p = 3.800 + 4.644/A 04 . The nonlinear coupling 
strengths of the a meson are given by 52 = —0.328 — 27.879/^4 0,4 (fm _1 ), and 
#3 = 38.862 — 184.191/A - 4 , whereas the self-coupling of the u field has the strength 
34 = 151.590 — 378.004/A ' 4 . For the pairing interaction, we take the strength of the 
delta function interaction as Vq = 343.7 MeV fm 3 , which was obtained by requir- 
ing that the experimental value of the proton pairing gap in 90 Zr (1.714 MeV) be 
reproduced with a given energy cutoff {E m3iX — A < 8.0 MeV). 

Because in many nuclei we have several solutions at different equilibrium defor- 
mations with similar energies, it is difficult to select the ground-state configuration 
uniquely. The procedure we employed is as follows. The basis deformation (3q is 
set equal to /?2m 5 following the results of Hirata et al., 4:! in which a constrained cal- 
culation^ was carried out for the quadrupole moment, Q201 to obtain the lowest 
minimum in the energy curve of each nucleus. 

The present calculation was performed by expansion in 14 oscillator shells for 
both the fermion fields and the boson fields. The convergence of this calculation has 
been tested with 20 shells for both the fermion fields and the boson fields. Following 
Ref.f 1 ^' we fix hw = AIA- 1 / 3 for fer mions. In what follows, we discuss the details 
of our calculations and the numerical results for the Zr isotopes. 
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Table I. The ground state properties of even Zr isotopes calculated with the parameter set TMA. 
Listed are the total binding energy, B to t, the binding energy per nucleon, B per , charge, neutron, 
proton, and matter root mean square radii, R c , R n , R p and R m , and the quadrupolc deformation 
parameter for the neutron, proton and matter distributions, /?2n, flip an d /?2m, with A the mass 
number and N the neutron number. 
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3.1. Binding energy per nucleon and two neutron separation energy 
The two neutron separation energy, S2 n , denned as 

S 2n (Z,N) = B(Z,N)-B(Z,N-2), (3-1) 

is quite a sensitive quantity to test a microscopic theory, where B(Z,N) is the 
binding energy of nuclei with proton number Z and neutron number N. The two 
neutron separation energy becomes negative when the nucleus becomes unstable with 
respect to two-neutron emission. Hence, the drip line nucleus for the corresponding 
isotope chain is the one with two less neutrons than the nucleus at which S^n first 
becomes negative. 
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Fig. 1. Binding energy per nucleon, B per , for even Zr isotopes as functions of mass number A 
obtained from the deformed RMF+BCS calculations (squares), the deformed RMF calculations 
(circles), the spherical RCHB calculations"*' (up triangles) and the experimental datcP^* (down 
triangles) . 
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Fig. 2. Two neutron separation energies, 5*2n, for even Zr isotopes as functions of mass number A 
obtained from the deformed RMF+BCS calculations (squares), the deformed RMF calculations 
(circles), the spherical RCHB calculations^" (up triangles) and the experimental data^ (down 
triangles) . 
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In Figs. Q and El we plot the results for the binding energy per nucleon of 
80-i08 2 r an( j resu lts for the two-neutron separation energies for the entire chain 
of Zr isotopes covering the proton and neutron drip lines. The figures also display the 
results of the RMF calculations, the results of the spherical RCHB calculations^) 
and the available experimental datal^U First, from Fig. ^ we can see that the 
RMF+BCS calculations give a better description of the binding energy per nucleon 
than the RMF calculations. The largest difference between the RMF+BCS results 
and the experimental values is less than 0.04 MeV. Noting here that most of the nuclei 
are deformed and that despite this fact, we did not readjust any parameter for our 
calculations, the agreement is quite remarkable. Second, from Fig. EJ we see basically 
good agreement between experiment and the present calculation. The values of 
S2n for the so-called giant halos^J ( 124 Zr— 138 Zr) are reproduced accurately. The 
strong variation in the experimental separation energy at the neutron magic number 
N = 50 is well accounted for by the present calculation. The small staggering for 
8 Zr, 88 Zr, 100 Zr and 114 Zr can be attributed to a mixture of the pairing interaction 
and the deformation effect. With the assumption of spherical shapes, it is predicted 
in Ref.' 2 -' that the drip-line nucleus for Zr isotopes is 140 Zr. We also obtain 140 Zr 
as the drip-line nucleus. In both Figs. ^ and [21 the deformation effect is clearly 
seen. The deformed calculations describe the experimental data much better than 
the spherical calculations. This once again points out the need for an appropriate 
relativistic calculation with both deformation and proper pairing interaction taken 
into account in order to obtain a reliable description of all the nuclei from the proton 
drip line to the neutron drip line. 

3.2. Root mean square neutron radii 

The root mean square neutron radius is another basic important physical quan- 
tity to describe neutron-rich nuclei. In the mean field theory, the root mean square 
(rms) neutron radii can be directly deduced from the neutron density distributions, 

Pn- 

'J p n r 2 dr~ 



Rn = K 



2\l/2 



/ Pndl 



1/2 

(3-2) 



In Fig. |SJ the root mean square neutron radii for Zr nuclei are presented. Al- 
though the calculation was done with the assumption of spherical shapes, the results 
of RCHE 20 1 are also shown for comparison. Two interesting features are clearly seen. 
First, the so-called giant halos ( 124 Zr— 138 Zr)22» are obtained. Second, nuclei with 
large absolute values /3 2m , 78 Zr- 82 Zr and 94 Zr- 112 Zr (see also Table. Hand Fig. H|, 
tend to have larger rms neutron radii. This can also be seen from the fact that some 
RMF results (for 100 Zr and 126 Zr) are larger than the RMF+BCS results. 

Another point to note here is that despite the second effect mentioned above, for 
the so-called giant halos, the present work gives relatively small values for the root 
mean square neutron radii. We believe that this is due to the harmonic oscillator 
basis we used in the calculations of the deformed nuclei. This deficiency will be 
resolved in the near future, although a great deal of effort is necessary to change 
the numerical method. Nevertheless, despite the small discrepancies for giant halos, 
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Fig. 3. The root mean square neutron radii for even Zr isotopes as functions of mass number A 
obtained from the deformed RMF+BCS calculations (squares), the deformed RMF calculations 
(circles), and the spherical RCHB calculation^^ (triangles). 



the present work provides good agreement with the RCHE 20 ' results and can give 
a reliable prediction of the rms neutron radii for all the nuclei from the proton drip 
line to the neutron drip line. 

3.3. Single-particle states and their occupation probabilities 

One interesting feature of exotic nuclei is the contribution from the continuum 
due to the pairing correlations. In this context, it is very interesting to study the 
amount by which the contributions from the continuum differ for the calculations 
with a constant pairing interaction and the calculations with a delta function inter- 
action. In Fig. |IJ we plot the occupation probabilities of 124 Zr for the neutron 
levels near the Fermi surface, i.e. in the interval —10 MeV < E s _ p _ < 8 MeV. We 
present the occupation probabilities of neutron single-particle states for two cases, 
the delta function interaction and the constant pairing interaction usually used in 
the BCS framework. In the constant pairing calculation, we used G n = 12.0/yl 
and G p = 30.0/^4, where A is the mass number of 124 Zr, and the pairing window 
£j — A < 2(41^4 -1 / 3 )P) Here, we used different pairing strengths, G p and G n , for 
protons and neutrons in order to obtain similar gap energies, A, than those obtained 
with the delta function interaction. We note that we use different pairing strengths 
for protons and neutrons because the level density around the Fermi surface is small 
for protons and large for neutrons, due to the continuum. The results obtained from 
the calculations using the delta function interaction and the constant pairing interac- 
tion are plotted in the left and right panels, respectively. Because in our calculation 
124 Zr is almost spherical, we denote the corresponding spherical quantum number in 
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Fig. 4. The occupation probability is represented by solid horizontal bar for each neutron single 
particle state of 124 Zr. The left panel displays the results of the RMF+BCS calculation with 
the delta function interaction, while the right panel displays the results of the constant pair- 
ing calculation. The dotted horizontal line represents the Fermi energy. In both panels, the 
occupation probabilities of the continuum states are multiplied by a factor of 5 for clarity. 
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Fig. 5. The wave function f(z,r±) for four typical resonant states in 124 Zr as functions of r± 
at z = 0.422. The corresponding energy (in MeV) and occupation probability are, respec- 
tively, for 13/2+ [606], E = 5.276 and v 2 = 0.005; for 9/2" [505], E = 1.857 and v 2 = 0.041; for 
7/2" [503], E = 0.464 and v 2 = 0.196 and for 5/2" [503], E = 1.683 and v 2 = 0.028. The Nils- 
son quantum numbers in square brackets are determined using the dominant component in the 
expansion of this wave function in terms of the anisotropic oscillator basis. 
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the left panel. We can see all the spherical shell model states with small spin-orbit 
splitting in the positive energy region. They are all resonance states supported by 
the centrifugal potential. A recently performed resonant RMF+BCS calculation^ 
yielded a quite similar spectrum for 124 Zr. It is thus seen that the contributions from 
the continuum states are very important. We see clearly that the BCS calculation 
with the delta function interaction can include more of the continuum states that are 
more localized inside the nuclear region, which correspond to resonance states. In 
addition, the delta function interaction includes a smaller contribution of the contin- 
uum states that are less localized in the pairing calculations. To see this point more 
clearly, we plot the corresponding resonant wave function / in Fig. [5J It is easily 
seen that these states do have similar behavior as bound states. On the other hand, 
in the constant pairing case, the occupation probabilities decrease monotonically as 
functions of the single-particle energy. 

We next point out an interesting feature of the vacancy between E ~ 2MeV 
and E ~ 5MeV in the single-particle spectra in the continuum seen in both calcu- 
lations. If we performed calculations for the single-particle states in the coordinate 
space using box boundary conditions instead of the harmonic oscillator expansion 
method, we would obtain many states in the continuum. These are the so-called 
scattering states, which have small probabilities in the nuclear region. We do not 
obtain these scattering states, at least in the region of several MeV excitation energy 
in the continuum. Hence, we find that the single-particle states near the continuum 
threshold do have large probabilities in the nuclear region and contribute to the 
pairing correlations in nuclei close to the neutron drip line. 

This unique feature of the delta function interaction in the BCS method is 
essential for the study of drip- line nuclei, where the Fermi energy is close to the 
threshold of the continuum. In this case, we have to estimate correctly the coupling 
between the bound states and the continuum states in order to pick up the resonance 
states, which have large amplitudes in the nuclear region. This would justify the use 
of such a simple RMF+BCS model to study all the nuclei, including the unstable 
ones from the proton drip line to the neutron drip line, as already demonstrated for 
the spherical case by Yadav et alP 

3.4. Effect of pairing on the deformation 

The quadrupole moments of proton, neutron and nucleon distributions are cal- 
culated as 

Qi = ^{^(6)), (3 . 3) 
= {2z 2 -x 2 - y 2 )i, 

where Y\ m represents the spherical harmonics, with I being the multi-polarity. The 
index i = p,n,m denotes the expectation value with respect to the proton, neutron 
and nucleon distributions, respectively. 

The deformation parameters are defined in terms of the liquid drop model with 
uniform density. We expand the sharp surface of the liquid drop as 

Ri(9) = R [l + /3 2i Y 20 (6) + (3^(9)} , (3-4) 

where R$ = 1.2A l l^ fm, in terms of the spherical harmonics under the assumption 
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of axial symmetry. The quadrupole moment of the nucleon distribution in the liquid 
model is calculated as 

Q% = ^^fom, (3-5) 

dropping terms of higher order in f3\ m . The deformation parameter fom is deter- 
mined such that Q l m reproduces Q m in the RMF calculation, and is given by 



a — V 5-7T 1 n 

P2m — — q T~nl 



?I^o ( 3 . 6 ) 

V57T 1 /o y 2 „2 ,,2\ 

"^~ARf (2Z " X _2/ )m - 



The deformation parameters fop and /?2n are similarly given by 



hp JL 3 — ^2^P' ( 3 ' 7 ) 



fon ~ ^J^Qn- (3-8) 

In the procedure to extract (3\i described above, we retain the linear relations 
between the moments and the deformation parameters so that it is easy to reproduce 
the moments from the deformation parameters in Table. [I] Regarding the determi- 
nation of the deformation parameters with terms of higher order in (3\i > we refer the 
reader to the description in RefP 




Mass Number A 



Fig. 6. The quadrupole deformation, /?2m, for even Zr isotopes obtained from the deformed 
RMF+BCS (squares) and the deformed RMF (circles) calculations as functions of mass number 
A. 
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We plot in Fig. El the quadrupole deformation, ftmj obtained from both the 
deformed RMF calculation and the deformed RMF+BCS calculation. It is easily 
seen that the pairing effect reduces the deformation of nuclei. More specifically, for 
some nuclei with \P2m\ < 0-2 as found from the RMF calculation, (i.e. 88 Zr and 
118 Zr— 124 Zr), the deformations are removed almost completely by the pairing effect. 
For other largely deformed nuclei (with ftm > 0.3), ft™ is reduced somewhat, but 
not as much as in the case of their weakly deformed counterparts. In Fig. we 
compared our predictions for the quadrupole deformation parameter fi\ with the 
empirical values.^ Except for 82 Zr, in which case our result is larger than the 
empirical value, the results are quantitively in good agreement with the empirical 
values. 
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Fig. 7. The square of the quadrupole deformation, /3| p , for even Zr isotopes obtained from the 
deformed RMF+BCS calculations (solid squares) and the empirical valuea"" (empty squares) 
as functions of mass number A. 



§4. Conclusion 



We have formulated the RMF theory with deformation and pairing. Conven- 
tionally, calculations for deformed nuclei are carried out using the expansion method 
in terms of the harmonic oscillator wave functions with pairing correlations treated 
using a constant pairing interaction with a pairing window. This method is, however, 
not applicable to the case of nuclei close to the neutron and proton drip lines, due to 
the importance of the resonance states in the continuum for such nuclei. In order to 
extend the RMF method so that it is applicable in these regions also, we introduced 
a delta function interaction for the pairing interaction, which allows the model to 
pick up resonant states by making the pairing matrix elements state dependent. The 
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delta function description has been demonstrated to work for spherical nuclei. 

In this paper, to demonstrate the applicability of the method, we have studied 
the Zr isotopes from the proton drip line to the neutron drip line. We calculated 
the binding energies, nuclear radii and deformation parameters of these nuclei. We 
also calculated the binding energy per nucleon and the two-neutron separation en- 
ergy. We found that the agreement with experimental values is very satisfactory. 
The effect of deformation is clearly seen in these quantities as the large variation 
of the deformation as a function of mass number. We found that the neutron radii 
increase monotonically with the neutron number, with some anomalies, where the 
deformation changes suddenly. In our results, the so-called giant halo effect is pre- 
served for nuclei close to the neutron drip line. The comparison of our results for 
02 P with the empirical values shows that our predictions for deformation parameters 
are quantitively in good agreement with the experimental results. 

The occupation probabilities in the continuum are important for the purpose of 
determining if the present pairing method is effective in the treatment of the pairing 
in the continuum. While the occupation probabilities decrease monotonically as the 
states deviate from the Fermi energy in the case of a constant pairing interaction, 
they exhibit characteristic behavior for the case of a delta function interaction. The 
occupation probabilities are large for those states whose wave functions have large 
overlap with the wave functions below the Fermi surface. It would be interesting 
to compare our results with those of the Hartree-Bogoliubov calculations for these 
nuclei. 

The deformation varies greatly for the Zr isotopes. The pairing correlations 
have the effect of reducing the deformation. This is clearly seen in our calculations. 
For those nuclei that have small deformations when calculated without pairing, the 
pairing correlations cause the nuclei to be spherical. 

In conclusion, we have carried out calculations for the Zr isotopes to demonstrate 
that the presently considered RMF method with deformation and pairing correlations 
is effective even for nuclei close to the drip lines. It would be very interesting to make 
similar calculations for many nuclei in all mass regions with the present method. 
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